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Abstract 

In this paper we examine the properties of U( 1) gauged Q-balls in two models with 
different scalar field potentials. The obtained results demonstrate that in the general case 
U{ 1) gauged Q-balls possess properties, which differ considerably from those of Q-balls 
in the nongauged case with the same forms of the scalar field potential. In particular, 
it is shown that in some cases the charge of 17(1) gauged Q-ball can be bounded from 
above, whereas it is not so for the corresponding nongauged Q-ball. Our conclusions are 
supported both by analytical considerations and numerical calculations. 


1 Introduction 

A simplest generalization of nontopological solitons, initially proposed in [1] and known as 
Q-balls [2], from the global £7(1) symmetry to the gauge 17(1) symmetry was proposed and 
analyzed in the pioneering paper [3]. Later this subject was examined in the well-known paper 
|4], in which gauged Q-balls (for simplicity, from here on, we call £7(1) gauged Q-balls “gauged 
Q-balls”) were examined analytically and numerically. One can also recall papers mum, 
where gauged Q-balls were examined mainly from a theoretical point of view, as well as papers 
iMEnunjiia, where solutions for gauged Q-ball were obtained numerically. 

It is clear that if the backreaction of the gauge field on the scalar field is small, the character¬ 
istics of gauged Q-balls (charge, energy, etc.) do not differ considerably from those of ordinary 
nongauged Q-balls (see paper [7], in which this issue was examined in detail). Meanwhile, the 
most interesting cases are those in which backreaction of the gauge field cannot be neglected. 
In the general case, the latter realizes not only if the corresponding coupling constant is large, 
but even when the coupling constant is rather small, but the other parameters of the solution 
are such that contribution of the gauge field is considerable. In the present paper we examine, 
both analytically and numerically, such differences between nongauged and gauged cases. In 
particular, we show that in some cases gauged Q-balls may exist for such values of the param¬ 
eters, for which nongauged Q-balls do not exist at all. And vice versa, it is possible that there 
are no gauged Q-balls for the values of the parameters for which nongauged Q-balls exist. The 
obtained results demonstrate that gauged Q-balls possess properties, which can be completely 
different from those of Q-balls in the nongauged case. 

The paper is organized as follows. In Section 2 we present the general setup and introduce 
the notations that will be used throughout the paper. In Section 3 we present some analytical 


1 


results related to gauged Q-balls. In Section 4 we examine numerically two models with different 
scalar field potentials. The obtained results are briefly discussed in the last section. 


2 Setup 

We consider the action, describing the simplest U( 1) gauge invariant scalar field theory in 
four-dimensional space-time, in the form 

S = I d 4 x - ieA + ieA,,cp) - V ( 0 * 0 ) - \f, v F^ ( 1 ) 


and take the standard spherically symmetric ansatz for the scalar and gauge fields P®: 

df(r) 




f{r) Ir^oo -t 0, 

A 0 (t,x ) 

= M r ), 

A)( r ) r->oo ~t 0, 

Ai(t,x ) 

= 0 , 


'W and /(r), 

A 0 (r) are real functions. Below 


dr 
dA 0 (r) 


= 0, 


r=0 


dr 


= 0, 


r=0 


( 2 ) 

(3) 

(4) 


the function /(r) has no nodes. Without loss of generality we take /(0) > 0. 
It is obvious that, according to ([2])-(j4]), we can use the effective action 


S e{i = 4t r / r~dr ( (w + eA 0 ) 2 f 2 - d r fd r f -V(f) + -d r A 0 d r A 0 


(5) 


where V(f) = V(<f)*(j>). For the scalar field potential, the conditions 


1/(0) = 0, 


dV 

df 


= 0 


( 6 ) 


f-o 


are supposed to fulfill in order to ensure the existence of the vacuum solution f(r) = 0, A 0 (r) = 
0. The equations of motion, following from effective action (]5j), take the form 


2 e(u/ + eA 0 )f 2 = --^( rA 0 ), 
r dr z 

. 9 „ Id 2 . 1 dV 

(W + eAo) f+ rd^ {rf) ~ 2? = 0 ' 


(7) 

( 8 ) 


For the numerical analysis, it is more convenient to use the combination a(r) = u + eAo(r) 
instead of the field A 0 (r). With this notation, equations dTJ), (jHJ) can be rewritten as 


.,2 f 2 


1 d 2 

« 2 / + --r^( r f 

r dr z 

where ^lki| _ = 0. The value of the frequency u is now defined as cu = lim a(r). 


Id 2 , , 


rdr^ ra ’ 

(9) 

ldV 


-= 0 

2 df ’ 

(10) 
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The charge of a gauged Q-ball can be defined aqj 


Q = 8n J (oj + eA 0 )f 1 2 r 2 dr = 87 T J a f 
0 0 


2 - 2 dr. 


( 11 ) 


According to El S3, the sign of a = u + eA 0 always coincides with the sign of oj, whereas A 0 = 0 
for oj = 0. Thus, without loss of generality we can consider oj > 0, which leads to Q > 0. The 
energy of a gauged Q-ball at rest is defined by 


E = in J [a 2 f 2 + d r fd r f + V(f) + d r ad r a ) r 2 dr. 
0 


( 12 ) 


It is well known that for ordinary (nongauged) Q-balls the relation ^ = oj holds. In [7] 
it was shown that the same relation also holds for U( 1) gauged Q-balls. We will use it for an 
extra check of our numerical results. 


3 Analytical considerations 


To begin with, let us discuss the allowed values of the frequency oj. Suppose that our scalar 
field potential is such that the relation 


]_dV 


= M 2 


f =0 


(13) 


holds. In the most cases nongauged (i.e., with global 17(1) symmetry) Q-balls exist only for 
oj < M, whereas Q —* 00 for oj — * M. So, it is universally accepted that for gauged Q-balls the 
values of the frequency oj are also bounded from above as oj < M (recall that we take oj > 0 ). 
Indeed, it was shown in [3] using the perturbation method in the effective coupling constant 
that the total energy of gauged Q-ball diverges for u = M even in the special case in which the 
corresponding nongauged Q-ball exists and has finite charge and energ}@. In [1] it was stated 
that the condition u < M is required to have localized solutions without oscillations for the 
scalar field. However, our numerical analysis shows that gauged Q-balls with finite charge and 
energy for oj — M may exist even if the corresponding solution in the nongauged case does not 
exist at all. Below we will present some analytical considerations which support this statement. 

First, we consider the usual case oj < M and suppose that there exists a gauged Q-ball 
solution with finite charge and energy such that A 0 (r) —> — ^ for large r. For such large r the 
equation for the scalar field can be rewritten as 


(cu 2 — M 2 )f — 2ue Q f + ~~r~ 9 ( r /) ( 14 ) 

47 rr r dr- 

1 The physical charge is defined by Q p h ys = eQ , but below we will use the charge Q defined by ( fill ) , which 
simplifies comparison with the nongauged case. 

2 Note that the correction to the background nongauged solution in [3] grows with r, which indicates the 
breakdown of the linear approximation at some r (moreover, as it was shown in [7], such a breakdown of the 
linear approximation for the correction is inherent to models of gauged Q-balls). Thus, the divergence of total 
energy in the linear approximation cannot be used as an indication of the absence of a solution to the full set 
of nonlinear equations. 
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The solution to this equation, tending to zero as r —> oo, can be easily obtained and takes the 
form 


/(r) = Ce 




u>e 2 Q 


47 Ty/W 2 


co- 


:, 2 , 2 , 


(15) 


where C is a constant and U(b,c,z) is the confluent hypergeometric function of the second 
kind, ft is not difficult to show (see Appendix A) that for \JM 2 — ui 2 r > 1 we get for (fT5|) 


f(r) ~ 


e -VM‘ 2 -ui 2 r 


1 + 

r 


uie 2 Q 
4 7T \/ M “ — 


(16) 


. . . —\/iP — r i i . rr • i 

This formula resembles the naively expected result /(r) ~ ----, the difference is caused 

by taking into account the electromagnetic potential Aq{t) —> — ^ for large r. But formula 

m has a singular behavior in the limit u — > M, which clearly indicates that the case u = M 

should be considered separately. 

Equation (fT4l) for u — M also has a solution, ensuring the finiteness of charge and energy. 
Indeed, let us suppose that solution to the full set of nonlinear equations exists and has the 
charge Q. Then, far away from the center of the Q-ball, we can write for the scalar field 


2Me 2 Q \cP 
~T^~ f+ rd^ <rf] 


0 . 


The solution to this equation, tending to zero as r —y oo, takes the form 


(17) 


K i 

/ (r) = C— 



(18) 


where C is a constant and K\ (b, z ) is the modified Bessel function of the second kind. For large 

r this solution behaves as _ 

J imJq ' 

e v 7i" 

fir) - 3 -■ (19) 

f*4 

On sees that due to the long-range action of the electromagnetic potential Ao(r), the behavior 
of the scalar field at r —> oo differs considerably from the nongauged case, in which one expects 
fir) ~ 1 for uj — M. Moreover, one may naively expect that the repulsive nature of the 
electrostatic interaction would prevent from forming a gauged Q-ball for such a value of u> 
(recall that usually the corresponding nongauged Q-ball has an infinite charge for oj = M). 
The argumentation presented above shows that it is not so. This statement is also confirmed 
by the numerical results, which will be presented in the next section. 

ft should be noted that it is possible to show analytically that in the limit oj —> M solution 
(fT5j) transforms into solution (118ft . see Appendix B. 

A few words about the case uj > M. The form of equation (fT4j) suggests that the corre¬ 
sponding solutions for the scalar field are oscillatory for r —y oo, leading to infinite charge and 
energy. More precisely, the leading term of a solution to equation (TT4|) takes the form 

C\ - cos (Vu) 2 — M 2 r- ^ = In (Vco 2 — M 2 r 

r V 4v ry/uJ 2 - M 2 V 

+C -2 - sin ( y/co 2 — M 2 r - -== 3 = In ( \/uj 2 — M 2 r 

r V 4ttVw 2 - M 2 V 
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The latter formula supports the assumption that there are no gauged Q-balls for u > M, which 
was also confirmed by the numerical analysis. 

One can think that it looks rather strange that there exist solutions with finite charge for 
oj < M, whereas no solutions with finite charge are expected for oj > M. However, a similar 
situation can be observed in the nongauged case. Indeed, let us take a scalar field potential of 
the form 

V{f) = M 2 f 2 -\\f\ N . (20) 

It was noted in [I3j that Q — > const 7 ^ 0 in the limit u — > M for N = ^ (one can easily 
check numerically that the corresponding solutions for the scalar field indeed exist). More 
precisely, Q ~ u for oj < M. Meanwhile, it is possible to show analytically that there are no 
Q-ball solutions for a j > M (to show it in a simple way one can use the scale transformation 
technique, proposed in JIT], supplemented by transformations of the fields themselves; see, for 
example, |TS], where such a method was applied to the more complicated case of gauged Q- 
balls). Thus, this example indicates that the situation with the existence of a finite charge for 
uj — M in the gauged case is not so unique. 

4 Explicit examples of U( 1) gauged Q-balls 

The numerical solutions for gauged Q-balls, which will be presented below, were obtained in 
two steps. On the first step the shooting method, which solves the boundary value problem by 
reducing it to the solution of the initial value problem, was used. According to this method, 
one should adjust the initial data for the system of equations at one of the boundaries (at the 
origin of gauged Q-ball in our case) in such a way that the solution satisfies required conditions 
at the second boundary (for large r in our case). The shooting method is very simple and 
it is easy to implement it, but it fails to find gauged Q-ball solutions for large r, where the 
scalar field falls off to zero exponentially, see relations HIB]) . (I19p . To find a solution for such 
large values of the coordinate r using the shooting method, one has to fine tune the initial 
data with very high accuracy, which may even exceed the truncation error of double precision 
floating-point numbers. Moreover, by taking the box of a small size one can mistake spurious 
solutions (such as oscillating solutions for u> > M ) for correct monotonic solutions of the 
boundary value problem. So, these problems make the shooting method not fully applicable 
for our task. In order to overcome them, the solutions obtained on the first step (i.e., for the 
values of r which are not very large) were supplemented by analytical solutions defined by 
formulas (fTBj) . (THU) . On the second step, the resulted “combined solutions” were used as the 
initial approximation of the solutions to the boundary value problem for discretized version of 
the system of equations ®,m- Then these solutions were improved iteratively by the Gauss- 
Seidel red-black relaxations, accelerated with the help of multigrid technique. The additional 
cross-checks of the final results were performed using the known theoretical relations for gauged 
Q-balls, such as, for example, ^ = (yy-ry) = uj [7]. 
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4.1 Model with 0 4 scalar field potential 

At first, we consider the potential of form 

V(f) = M 2 f 2 - A/ 4 , (21) 

where A > 0. At the very beginning it is convenient to make the following redefinition of the 
coordinate r and the fields: 

R=Mr, G(S)A«M, F(R) = ^f(r). (22) 

In these notations, the system of equations ([9]) , (TTOl) takes the form 

2a 1 GF 2 = 4 (RG)", (23) 

R 

G 2 F + ^-(RF)" — F + 2F 3 = 0, (24) 

R 

where ' = a\ = y. We see that the only effective parameter in this system of equations 
is oq. Since M 2 , w e will be looking for solutions such that G(oo) < 1, which 

corresponds to oj < M. 




tn/M oo/M 


Figure 1: Q(u) for different values of the parameter a± (thick lines). The thin lines stand for 
the nongauged case. The circles on the plots mark the points with — 1. 

The charge of the Q-ball takes the form 

OO 

Q = j 8vr J GF 2 R 2 dR = jQ, (25) 

o 
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whereas the energy is 

00 

E = ^- 4tt [ (G 2 F 2 + d R F8 R F + F 2 - F 3 4 + —d R Gd R G] R 2 clR = ^-E. (26) 

A J V 2«i / A 

0 

In Fig. Q] one can see several examples of Q(u) diagrams (expressed in the dimensional 
variables Q and |4). Since in the nongauged case there exists a solution for oj = 0 with the zero 
charge [IS], the gauged solution for u = 0 simply coincides with it [3], which explains why the 
curves start from the point oj = 0, Q = 0. Meanwhile, in the nongauged case the Q-ball charge 
tends to infinity while —* 1 [16]. As it was demonstrated in Section 3, it is not so for the 
gauged case, in which one may expect the existence of a Q-ball with finite charge and energy. 



Figure 2: ^j(Q) for different values of the parameter a.\ (thick lines). The thin lines stand for 
the nongauged case. The circles on the plots mark the points with j? = 1. 


It is confirmed by the plots in Fig. |T] demonstrating the existence of Q-balls with finite charges 
for — 1 for different values of the parameter a^o We also performed the numerical analysis 
for jj > 1. No solutions of form ©-(HD with finite charge and energy were found. 

In Fig. [2] one can find the ^(Q) dependencies for different values of the parameter a\ in 
comparison with the nongauged case. Profiles of the scalar held for different values of R and 

3 The numerical analysis was also performed for the values of the parameter a\ larger than those used in 

Fig- CD However, no essential changes in the behavior of the Q(u>) dependencies were found. 
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Figure 3: Profiles of the scalar field for different Figure 4: Comparison of the scalar field pro¬ 
values of Here ^JoT\ = 0.05. files (thick lines) with the asymptotes defined 

by Eqs. (fT6j) . (TT9]i (thin lines) for different val¬ 
ues of jj. Here yfa[ = 0.05. 


for yfal = 0.05 are presented in Fig. [3J We have also performed a comparison of the numerical 
solutions for the scalar field at large R with the asymptotes defined by Eqs. (fTfilh (fTOlh The 
result is presented in Fig.|4l demonstrating a remarkable agreement with theoretical predictions. 


4.2 Model with the piecewise parabolic scalar field potential 

Now we consider the piecewise potential of form 

V(f) = M 2 f 2 0^1 - + M 2 v 2 e - lj , (27) 

where 0 is the Heaviside step function with the convention 0(0) = | (such piecewise potentials 
for the case of nongauged Q-balls were introduced in pQ and thoroughly examined in [17, T8]). 
Again, it is convenient to pass to the new variables 

R = Mr, G(R) = jja(r), F(R) = ±f(r). (28) 

In these notations, the system of equations (jHD, (flUl) takes the form 

2 a 2 GF 2 = hRG)", (29) 

it 

G 2 F + ^{RF)" - F6 (1-F 2 ) =0, (30) 

2 2 

where «2 = which is the only effective parameter in this system of equations. Again, since 
jf~j/ r \f =0 = Af 2 , we will be looking for solutions such that G(oo) < 1, which corresponds to 
co < M. 

The charge of the Q-ball takes the form 

OO 

Q = W 81 f GF2,fdR = Jp <?• ( 31 ) 

0 


8 


















whereas the energy is 


OO 

E= Jj An J (° 2F2 + d * Fd * F + F 2 9(l-F 2 ) + 9 (F 2 - 1) + J-d R Gd R G) R 2 dR (32) 


= — E. 
M 


In Fig. [5]one can see several examples of Q(u) diagrams for this model. We see a completely 




0 0.2 0.4 0.6 0.8 

co/M 


0 0.2 0.4 0.6 0.8 1 

co/M 


Figure 5: Q(u) for different values of the parameter (thick lines). The thin lines stand for 
the nongauged case. The circles on the plots mark the points with = 1, the triangles mark 
the points with ^ = 0, the asterisks mark the points with ^ = oo, the boxes also mark the 
points with = 1. 

unexpected behavior of the corresponding Q(u) dependencies. First, we see that, contrary to 
the case of the previous model, now the parameter uj does not uniquely define the charge of 
the Q-ball. Indeed, except the Q-ball corresponding to a minimal value of u (these points are 
marked by asterisks), there exist two Q-balls with different charges for each value of uj. We 
have the following explanation of this fact. 

Indeed, in the nongauged case the charge tends to infinity in the limits uj — » 0 and uj —>- M 
PSj. According to the results of Section 3, the nongauged limit uj —> M transforms into the 
Q-ball with uo = M. This Q-ball corresponds to the lower points with = 1 in Fig. 0 (these 
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points are marked by the circles), and its existence is not surprising. Now let us consider 
another limit, namely u — > 0. From [3] we know that if to = 0 for a gauged Q-ball, then A 0 = 0. 
Since there is not a Q-ball solution with the zero charge for oj — 0 in the nongauged case [18] , 
there should be no such solution in the gauged case too. It is also improbable that the charge 
of a gauged Q-ball tends to infinity in the limit ui —* 0 — the value of u tends to zero, whereas 
A 0 (r) < 0 is a monotonically growing function such that oo + eA 0 (r) > 0 for any r. Thus, 
|eA 0 (r)| < oj —* 0, whereas a solution for A 0 (r) should support the existence of a large charge. 
The latter situation seems to be unrealizable. So, the Q(u) curve modifies with respect to 
the nongauged case in order to overcome this problem and to maintain its continuity, which is 


realized in Fig. [5]— the curve just turns back 
At this point, the value of the charge behaves 



Q 


at the point of the minimal possible value of u. 
as Q* ~ ^ ~ for small values of ck 2 . 



Q 


Figure 6: ^(Q) for different values of the parameter a 2 (thick lines). The thin lines stand for 
the nongauged case. The circles on the plots mark the points with = 1, the triangles mark 
the points with ^ = 0, the asterisks mark the points with = oo, the boxes also mark the 
points with -^ = 1. 

It is interesting to note that now it is the parameter G(0) that uniquely characterizes a 
gauged Q-ball, not uj like in the nongauged case or even in the gauged case discussed in the 
previous subsection. We can also see that the curves in Fig. [5] become smaller while increasing 
the value of the parameter a 2 . For a 2 > 0.013 no gauged Q-balls were found. Analogous 
observation of nonexistence of gauged Q-ball for the values of the coupling constant larger than 
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some critical value was made in H2j, where the model with supersymmetry motivated scalar 
held potential, which has the form similar to (12711 . was examined. 

A remark is in order here. The results presented in this and in the previous subsections 
suggest that there exists a maximal possible charge of gauged Q-ball (of course, its value should 
depend on the model at hand). We think that it is not so in the general case; see, for example, 
till . However, this statement seems to be valid for models with oo. It should be 

L J by I / —^0 ' 

noted that this restriction is not connected with the restriction on the charge of stable gauged 
Q-balls [3], which was shown to be incorrect in [7j. 



0 10 20 30 40 50 

R 


Figure 7: Profiles of the scalar field for different values of jj. Here ydTj = 0.02. 


In Fig. |6]one sees the M (Q) dependencies for different values of the parameter «2 in com¬ 
parison with the nongauged case. The cusps on the corresponding curves, which are marked by 
triangles, correspond to the points with (p = 0, which are also marked by triangles in Fig. 0 
The existence of the cusps follows from the fact that (p = u for gauged Q-balls [7j. 

For the completeness, we also present the scalar field profiles (Fig. [7]), the effective charge 
density q(R) = GF 2 (Fig. |8]) and the effective pressure p(R) (Fig. [9]) for different values of 
p and for yTT^ = 0.02. The effective (dimensionless) pressure is defined in the standard way 
through the energy-momentum tensor [ 15 ] as 


p = G 2 F 2 


1 (dF\ 2 1 UG\ 2 

3 \dR) + 6a 2 \dR ) 


(F 2 9(l-F 2 ) +e(F 2 - 1)) , 


(33) 
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0 10 20 30 40 50 

R 



Figure 8: Profiles of the effective charge den- Figure 9: Profiles of the effective pressure 
sity q(R), = 0.02. P{R), \/®2 = 0.02. 


one can show that the equality J 0 °° p(R)R 2 dR = 0 (the so-called “von Laue condition” [TUj 1 
fulhlls in the general case (this equality was also used for the additional cross-check of our 
numerical results). The values of fj for these plots are chosen such that they approximately 
correspond to the marked points in Fig. [5] in the clockwise order: = 1 — circle, ~ 0.96 

— triangle, ps 0.46 — asterisk (the point of the turnover), fj = 1 — box. 

There exists a dip in the vicinity of R = 0 on the charge density curve in Fig. [8] for = 1, 
clearly indicating the repulsive nature of the electrostatic interaction in the gauged Q-ball. 
Such dips on the charge density curves appear for Q-balls from the upper part of the Q(u) 
diagram (Fig. [5]) for fj > 0.506. It is also interesting to note that, according to Fig. [9l the 
pressure can be negative even in the center of gauged Q-balls. 


5 Conclusion 

As it was demonstrated above, gauged Q-balls posses surprising properties which differ consid¬ 
erably from those of Q-balls in the nongauged case. Namely, it was shown that there may exist 
gauged Q-ball solutions with a finite charge even for uj = M (if ^y^| / ._ 0 = M 2 ^ 0), which 
is usually impossible in the nongauged case. The corresponding analytical considerations were 
supported by the numerical calculations in two models with different scalar held potentials. 
Moreover, numerical analysis shows that even when Q —> oo as uj —> 0 in the nongauged case, 
the charge in the gauged case remains finite for all allowed values of uj. As a consequence, the 
regions of allowed frequencies appear to be different from those in the nongauged case. 
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Appendix A: The behavior of f(r) for \JM 2 — oj 2 r 1 

Let us take the integral representation of the confluent hypergeometric function of the second 
kind U(b,c,z ), which has the form 


U(b,c,z ) = 


m 


e~ zt t b - 1 (l + ty-h-'dt. 


(34) 


For solution (TT5j) we get 


f{r) ~ e 


—y/M 2 —u) 2 r I —2V M 2 —u) 2 rt 


ujeQ 


cue z Q 


t 4?r V M 2 -uj 2 + t) 4ttV M 2 -uj 2 df 


(35) 


Let us change the variable t such that \JM 2 — oj 2 rt = t. We get 


f{r) 


\JM 2 —uj 2 r 


-2t 


y/M 2 — co 2 r + t 


47 r\/^2 — 


dt. 


(36) 


It is clear that, due to the exponential suppression, the main contribution to the integral in 
is achieved in the region t ~ 1. Thus, for y/M 2 — oj 2 r 1, we can rewrite the integral in 
as 


- 2 1 


t 


y/M 2 — u 2 r + t 


4? rV^ — ufi 1 


dt 


-2 1 


t 


yj M 2 — uj 2 r 


4 7 ry/~M^ — oj 2 


dt 


y/M 2 — LO 2 r 


OO 


4 7T V — cj 2 


dt. 


(37) 


Finally, for (l36j) we obtain (up to a constant depending on u>) 

p—y/M 2 —oj 2 r 

f(r) 


ue 2 Q 


T 4ttV M 2 — tu^ 1 


(38) 


which coincides with ffTfil) . 


Appendix B: Correspondence between the scalar field so¬ 
lutions for oj < M and oj = M 


Let us again take the integral representation of the confluent hypergeometric function of the 
second kind U(b , c, z). We will be interested in the limit uj —> M for a fixed r. In this case (1361) 
can be rewritten as 


fir) ~ - lim 

r lj->m 


-2 1 


y/ M 2 — uj 2 r + t 


uje z Q 


4 7r y/ — cu^ 1 


dt. 


(39) 
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In the limit u —> M the term 


V M 2 —oj 2 r-\-t 


can be represented as 


- 1 - 


yj M 2 — uj 2 r 


yjM 2 — uj 2 r + t yjM 2 — u 2 r + t 


1 - 


Thus, we get 


- 2 1 


lim / e , __ 

J \ y/M 2 — oj 2 r + t 


ujerQ 


Airs/ -AT ^ — 


dt « lim / e 2 M 1 — 
J 
0 


\J M 2 — l o 2 r 


yjM 2 — (jj‘“^ T \ 4 tt \/ M 2 —ui 2 


(40) 


dt (41) 


= lim / e 


-2 1 


yj M 2 — oj 2 r\ \/M 2 -u ) 2 r 


dt — e 2t e 4 ,rt dt. 


Now let us define y = ■ / 2Me2< ? 
can be rewritten as 


rp rp - / 

1 i x — \/ Me 2 Q\ 


■ t. In these notations, the last integral in (14111 


- 2 1 - 


Me^Q r 


e e 47 rt at = 


- --- oo 

dt = J Me ^ r Je-^Ddx 


(42) 


Me 2 Q r 


87 r 


e 2 ( x+1 /dx+ / e 2 ( x+ *)dx 


,0 1 

Now let us make the redefinition x —> - in the last integral in (142]) . We get 


Me 2 Q r 


87T 


e 2 { x+ *)dx + / e 21 = 


,-S 




It is convenient to introduce the new variable 2w = a; + -, which leads, according to the limits 
of integration, to x — w — yjw 2 — 1. The integral in (143]) can be rewritten as 

’Me 2 Qr 




87r 


' ——— ' dw 


y/w 2 — 1 


(44) 


Me 2 Q r 


2tt 


y e wy Vw 2 — 1 dre, 


where we have performed integration by parts in the last step. Recalling the integral represen¬ 
tation of the modified Bessel function of the second kind K x (y), which has the form 


Ki(y) = J e “V «> 2 - ld«j, 


(45) 


combining formulas 


and taking into account the definition of y , we arrive at 


f(r) 


2 Me 2 Q 


(46) 


which obviously corresponds to (1451) . 
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